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Abstract Neutron scattering has been used to inves- 
tigate the magnetic correlations in the the spin ice mate- 
rial dysprosium titanate, Dy2Ti207. An isotopically en- 
riched sample was used to minimise neutron absorption. 
In zero field no magnetic order was observed down to 
50 mK but the magnetic diffuse scattering was in qual- 
itative agreement with that expected for the disordered 
low temperature state of dipolar spin ice. Application 
of a field of « 0.8T in the [100] direction led to long- 
range order. With the field applied in the [110] direc- 
tion a coexistence of long range ferromagnetic and short 
range antiferromagnetic order was observed. This is at- 
tributed to the pinning of only half the spins by the field. 
The hysteresis loops in both field orientations displayed 
unusual steps and plateaus. 
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1 Introduction 

Ho2Ti 2 C>7 and Dy 2 Ti 2 07 are regarded as good examples 
of spin-ice materials n],0,||. The magnetic rare earth 
ions occupy a pyrochlore lattice, a cubic array of corner- 
linked tetrahedra. In the near neighbour spin ice model 
the lattice is occupied by (111) Ising spins that are cou- 
pled ferromagnetically to nearest neighbours [Q . The re- 
sult is a disordered, macroscopically degenerate, ground 
state that is determined only by the rule that two spins 
must point into and two out of each tetrahedron. This 
"two in, two out" rule is analogous to the "ice rules" 
that control the proton arrangement in water ice, and 
so the near neighbour spin ice model maps exactly onto 
Pauling's model of the proton disorder in water ice j|, 
p|. The near neighbour spin ice model qualitatively de- 
scribes much of the behaviour of the real materials, but 
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for an accurate description the dipolar spin ice model 
has been developed [0] . In this model, the near neighbour 
ferromagnetic coupling between (111) spins is dipolar in 
origin and it has been shown that the long range part 
of the dipolar coupling maintains a spin-ice like disor- 
dered ground state even in the presence of an antiferro- 
magnetic near neighbour super-exchange gj] . The dipolar 
spin ice model has been shown to give an accurate de- 
scription of the zero field magnetic neutron scattering of 
Ho2Ti2 07 g in the spin ice regime (T < 2 K) as well 
as the specific heat of Ho2Ti207 and Dy 2 Ti207 through- 
out the low temperature range (50 mK - 20 K) |£,£§. The 
dipolar coupling in Ho2Ti207 and Dy2Ti2C>7 is of similar 
magnitude (D nn w 2.35K for both Dy and Ho) while the 
antiferromagnetic exchange is weak in both cases (-1.2K 
for Dy and -0.52K for Ho) §,§. The main difference 
between these two materials appears to arise from the 
nuclear spins. Ho2Ti207 has a single isotope ( 165 Ho) and 
strong hyperfine coupling that complicates the spin ice 
freezing process below « 0.8 K J3]. Dy2Ti2C>7, on the 
other hand, has much weaker hyperfine coupling that 
does not affect the dominant spin ice behaviour in the 
temperature range of interest. It might therefore be con- 
sidered the simpler spin ice material. However, the nat- 
ural mixture of Dy isotopes leads to significant neutron 
absorption (994 barn for 2200ms -1 neutrons). Thus, al- 
though neutron scattering should be the ideal micro- 
scopic probe of the spin correlations and dynamics of 
this system, the most significant experimental observa- 
tions so far are bulk measurements [J3|. To minimise the 
absorption problem we have prepared a single crystal of 
162 Dy2Ti2C>7. In this paper we present the initial results 
of our neutron scattering study of this sample, in both 
zero and applied magnetic field. 



2 Experimental 

The absorption cross section of natural Dysprosium is 
<j a = 994 barn, the large value being due mainly to 
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the presence of 164 Dy. Isotopically enriched Dy2C>3 with 
composition 156 Dy< 0.01%, 158 Dy< 0.01%, 160 Dy= 0.02%, 
161 Dy= 0.47%, 162 Dy= 96.8%, 163 Dy= 2.21% and 164 Dy 
= 0.5% was prepared by Goss Scientific Instruments Ltd. 
Enrichment in 162 Dy was considered optimum as it has 
a high abundance and the dominant impurity is 163 Dy 
which has a slightly lower absorption cross section. En- 
richment in 163 Dy would have lead to 164 Dy as the domi- 
nant impurity. Through enrichment, the absorption cross 
section was reduced by a factor of five to 208 ± 13 barn. 

A single crystal of Dy2Ti20y was prepared from the 
isotopically enriched Dy2C>3 and TiC>2 by the floating 
zone technique || . Neutron scattering was carried out at 
the ISIS facility on the indirect geometry spectrometer 
PRISMA. It was configured in the diffraction mode so 
that sixteen 3 He tube detectors are used simultaneously. 
Rotation of the crystal allows a rapid mapping of a large 
section of reciprocal space, making PRISMA ideal for 
observing magnetic diffuse scattering. 

In the first experiment the crystal was cooled by a 
3 He sorption refrigerator for measurements in zero ap- 
plied magnetic field. It was aligned with [110] vertical, 
such that the scattering plane contained (hhl) wavevec- 
tors. In the second experiment an Oxford instruments 
7T vertical field cryomagnet, with dilution refrigerator 
insert, was used. For the applied field measurements two 
field orientations were studied: [110] as above and [100] 
corresponding to a (Okl) scattering plane. The data were 
normalised to monitor counts and vanadium to remove 
the characteristic flux profile of the spallation source. 
Since the sample is still significantly absorbing a correc- 
tion for this was also applied. 

3 Results 

With [100] vertical, a map of reciprocal space was made 
at the base temperature (« 70 mK) in zero field. The 
diffuse scattering maxima observed agree with those pre- 
dicted in the (hhl) scattering plane for the dipolar spin 
ice model. For example there is a diffuse feature at 3, 0, 
in (hkO) and 0, 0, 3 in (hhl) M. On application of a field 
the diffuse scattering disappeared and was replaced by 
magnetic Bragg peaks at the Q = positions. A field of 
~ 0.7T was sufficient to saturate these peaks. 

This behaviour is readily explained as the applied 
field breaks the degeneracy of the six "two in, two out" 
spin configurations of the elementary tetrahedron. The 
pyrochlore lattice can be described as a face centred cu- 
bic lattice with a tetrahedral basis, and the degeneracy 
breaking means that every tetrahedron adopts the same 
"two in, two out" state with a net moment in the direc- 
tion of the applied field, [100]. This non-collinear ferro- 
magnetic structure allows the observed magnetic Bragg 
peaks at the Q = positions such as 2, 0, 0. Interest- 
ingly the experimental magnetization did not develop 
smoothly, but in a series of steps. Hysteresis was ob- 
served on cycling the field. 




Fig. 1 Diffuse scattering in the (hhl) plane. In zero field at 
270mK no magnetic Bragg peaks are observed. All resolution 
limited intense features are of nuclear origin. 
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Fig. 2 Scattering in the (hhl) plane with a field of 1.5T 
applied on [110] at ~ 60mK. Magnetic Bragg peaks have ap- 
peared at positions such as 2,0,0 and the diffuse features ob- 
served in zero field at positions such as 0,0,3 have sharpened 
into features elongated on [00/]. 



The [110] direction is a hard direction of magnetiza- 
tion M. In zero field, the magnetic scattering was ob- 
served to be diffuse, and characteristic of the disordered 
low temperature state of dipolar spin ice B , see Figure 
[y. Application of a small field (w 0.1 T) in this direction 
again caused the appearance of the Q = Bragg peaks; 
however, unlike the [100] direction, the diffuse scatter- 
ing features did not disappear. As the field was raised 
up to 1.5T the diffuse scattering sharpened around the 
Q = X positions such as 0,0, 1 without becoming res- 
olution limited. This is seen clearly in figure |2| Similar 
features have been observed in the neutron scattering of 
Ho2Ti2 07 Q,0], adding strength to the idea that the true 
ground state of these systems is a Q = X structure that 
is dynamically inhibited from being accessed on experi- 
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Fig. 3 Integrated intensity of the 0,0,2 magnetic Bragg peak 
and the 0,0,3 diffuse feature (scaledx3 for this plot) as the 
field is scanned at ~ 70mK. 



mental timescales Effl. Again the magnetization versus 
field curve was observed to have several sharp steps and 
plateaus (Figure ||). 

The formation of the Q — X structure is consistent 
with the spin ice rules p|. In this field orientation, as- 
suming perfect (111) spins, only two of the spins of the 
tetrahedral basis have a component along the field direc- 
tion. These form "in-out" [110] chains parallel to the field 
forcing the remaining two spins per tetrahedra into "in- 
out" [110] chains perpendicular to the field, as illustrated 
in Figure [|. The perpendicular chains are not coupled 
by the spin ice rules which, in the absence of any fur- 
ther neighbour coupling, would lead to two dimensional 
[110] Bragg sheets of scattering extended along [00/] in 
the scattering plane. The diffuse features in the experi- 
mental pattern are indeed extended along this direction, 
but the sharp build up of intensity around the Q = X 
points indicates a strong tendency to prefer Q = X short 
range ordering of the perpendicular rods (see Figure g) . 
Neglecting interference between the scattering from the 
two spin sets (perpendicular and parallel to the field) one 
arrives at the conclusion that the Q = Bragg scatter- 
ing arises from the parallel rods and the Q — X diffuse 
scattering arises from the perpendicular rods. However, 
it may be a crude approximation to separate the scat- 
tering in this way. 

In conclusion, we have obtained accurate neutron 
scattering data for Dy2Ti2C>7 that is in qualitative agree- 
ment with theoretical expectations for a spin ice mate- 
rial |l],|2|,[|,[|. It is of interest that, even in a relatively 
strong field along [110], the system remains only par- 
tially ordered. It is also noteworthy that the magnetic 
hysteresis loop shows several steps and plateaus. An un- 
derstanding of these effects awaits a detailed study of 
the static and dynamic properties of the dipolar spin ice 
model [p| in an applied magnetic field. 




Fig. 4 A fragment of the pyrochlore lattice illustrating the 
Q — X structure. The applied field is projected into the plane 
of the page and pins the unshaded spins into ferromagnetic 
chains. The shaded spins are not coupled to the field but 
are governed by the "ice rules" . The shaded spins may form 
two structures, either Q = with all shaded spin chains 
parallel or Q — X with shaded spin chains antiparallel. This 
degeneracy is not raised by the field. 
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